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^ ■ Abstract 

Impact of the cone algorithm parameters E cut , E seec i, R on the efficiency and char- 
acteristics of the reconstructed jets in p — p collisions at the energy y/s = 200 GeV 
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is studied. The PYTHIA Monte Carlo generator is used for event generation. The 
dependence of dijet production fraction on the parton transverse momentum p± at 
different algorithm parameters is analyzed. The dependence of reconstruction ef- 
ficiency of parton energy in dijet events and two leading jets in iV-jet {Nj et > 2) 
events on E cut} E see d, R is studied. Monte Carlo results are compared with predic- 
Q-i| tions made in the framework of ^-scaling and experimental data obtained at RHIC. 

The independence of the slope parameter f3 of the scaling function, ip{z) ~ z _/3 , 
on the algorithm parameters over the energy range E^ et = 25 — 60 GeV is found. 
The strong dependence of the invariant cross section and the slope parameter on 
^ the algorithm parameters with decreasing of E^ et for E^ et < 25 GeV is observed. 

1 Jet-finding algorithm and jet properties 

Perturbative Quantum Chromodynamics (QCD) predicts the production cross sections 
for parton-parton scattering in p — p collisions at high px- The outgoing partons from 
the parton-parton scattering hadronize to form jets of particles. Calculations of high-pr 
jet production involve the folding of parton scattering cross sections with experimentally 
determined parton distribution functions (PDFs). Measurements of the inclusive jet cross 
section, the dijet angular distribution, and the dijet mass spectrum, are usually used to 
test predictions of perturbative QCD. 

We used the Monte Carlo code PYTHIA 5.7 for event generation. The mechanism of 
jet production involves the hard scattering process, the initial and the final state radi- 
ation and the multiple parton interaction pp. The dependence of the jet reconstruction 
efficiency, the parton transfers momentum P^ art accuracy and the inclusive cross sections 
of jet production on the algorithm parameters are studied. The z-scaling predictions for 
transverse jet spectra are compared with MC-results and RHIC data. 

1.1 Jet-finding algorithm 

We used the adapted by DO experiment algorithm for jet reconstruction [2]. The jet- 
finding algorithm consists of the following steps: 
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1. The particles with energy E l T > E seed are assigned as "seeds". 

2. Start with highest-Py "seed" and iterate over all "seeds". The "seed" -direction gives 
the first approximation of the jet direction (r] Jet , cf)j et ). 

3. All particles with distance Ri to the jet axis in the (rj, <p) space lesser than R are 
included into jet 

R l = y/(VJet ~ Vi) 2 + {<Pjet - 0,) 2 < R, (1) 

where (rji,^) are the particle direction. 

4. The energy E^ et and the direction of jet are calculated using ([2]): 
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5. The steps (3)-(4) are iterated until the jet-direction is stable. 

6. If Ej, et < E cut , the jet is discarded. 

7. Overlapping jets are merged or splitted depending on the energy fraction in the over- 
lapping region. 

8. Repeat these steps for all "seeds". 

9. The jet directions are recalculated using an alternative definition as given in (j3J) 







Jet 



tan 



(E^) 2 + (£*4) 



>Jet 



tan 



TE i 
TE i 



7] j et = -ln[tan(0 j et /2)] 
(3) 



where i corresponds to all particles whose centers are within the jet radius R, E % x = 
Ei sin Oi cos 4>i,E % y = Ei sin ^ sin 4>i,E\ = Ei cos 0j. 

In the paper we consider E see d, R, E cut as the parameters of the jet-finding algorithm 
and analyze their influence on characteristics of reconstructed jets. 



1.2 Jet reconstruction efficiency 

We define the jet reconstruction efficiency as a probability to find two or three jets 
in event. Figure [1] shows the probability to find two or three jets as a function of the 
transfers momentum p± of the hard process for different sets of parameters. The left 
frame is for the fixed E seec i, R and for different values of E cut , and the central frame is for 
the fixed E see d, E cut and for different values of R. One can see that values of the algorithm 
parameters E cut , R define the low limit P^ m of the pj_ spectrum and the ground boundary 
P^ at of saturation region for dijet events. For example, the dijet reconstruction probability 
drops very fast for p± < 13 GeV at E cut =7 GeV. The low limit of the p± spectrum 
increases with E cut and as R decreases. The probability to find two jets for given p± has 
a maximum at E cut w p±/2. It diminishes with decreasing of E cut . Part of events are 
reconstructed as three-jets. The right frame is for fixed E cut and R and different values of 
E see a. As seen from Figure[T]the influence of E see d on jet reconstruction is very small while 
E S eed « P±- Therefore, following results are presented for E see d=l GeV. The parameters 
E cut and R have an influence on the number of reconstructed jets. Therefore the following 
results are presented for dijet events and for two leading jets of iV-jet {Nj e t > 1). 

1.3 Parton transverse momentum reconstruction 

We characterize the parton transverse momentum P^ art reconstruction efficiency via 
the dependence of the mean jet transverse energy < E^ et > and the E^ et distribution 
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Figure 1: The dependence of the probability to find two or three jets in p — p collisions on 
the transfers energy of the hard process p± at the different values of the algorithm param- 
eters: (a) E seed =1.0 GeV, R=0.7, E cut =5, 7, 14 GeV; (b) E cut =7.0 GeV, E seed =1.0 GeV, 
R=0A, 0.7, 1.1; (c) E cut =7.0 GeV, i?=0.7, E seed =0.5, 1.0, 1.5 GeV. 



width (RMS) on parton transverse momentum P^ art . These dependencies are obtained 
from analysis of the dijet distribution events as a function of jet energy E^ et for narrow 
parton momentum bins. One of them is presented in Figure [2]^a). Figure [2](b) shows the 
dependence of the mean jet transverse energy as a function on P^ art for a fixed value of 
E cut and E seed and different values of R. The linear dependence of < E^ et > on parton 

Similar results are found for the 
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another sets of the jet finding algorithm parameters. 
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Figure 2: (a) The dependence of number of dijet events on the jet transverse energy E^ et 
for 20 < P^ art < 21 (GeV/c). (b) The dependence of the mean jet transverse energy 
< E^ et > on the parton transverse energy P^ art for two leading jets in iV-jets events 
(N Jet > 2) at E cut =7 GeV and R = 0.4, 0.7, 1.1. 



Figure [3] shows the dependence of RMS as a function of P^ art for fixed values of R 
and E seed and different values of E cut . The left frame is for dijet events. The slight growth 
of width is observed for P^ art > P^ at . We would like to note, that the momentum P^ at 
defines the low boundary of the saturation region of the probability to find dijet events. 
The diminution of the dijet reconstruction probability for P^ art < P^ at leads to decrease 
of RMS values. The lower E cut provides the smaller E^ et distribution width. The dijet 
reconstruction probability decreases in this case. The right frame is for two leading jet of 
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Figure 3: The dependence of RMS for dijet events (a) and two leading jets in iV-jets 
events (Nj et > 2) (b) on the parton transverse momentum P^ art at R — 0.7 and E cut = 
14, 7, 5 GeV. 



iV-jet events, where iVj et > 1. The jet transverse energy distribution width is proportional 
to P£ art in this case. The parton transverse momentum reconstruction efficiency is almost 
the same for different values of E cut . 




Figure 4: The dependence of RMS for dijet events (a) and two leading jets in iV-jets 
events (Nj et > 2)(b) on the parton transverse momentum P^ art at E cut =7 GeV and 
R = 0.4,0.7,1.1. 

Figure H] demonstrates the dependence of RMS on P^ art for fixed values of E cut and 
E see< i and different values of R. The left frame is for dijet events. The growth of width for 
P^ art > P^ at and alteration of RMS for P^ aTt < Pj, at are found. The higher R provides 
the smaller E^ et distribution width. The difference of RMS for R = 0.7 and R = 1.1 is 
found to be not so large. Results of analysis with higher R are potentially more sensitive 
to background. The value of the RMS raises substantially for R=0A. The right frame 
is for two leading jet of iV-jet events (Nj et > 1). The jet transverse energy distribution 
width is proportional to P^ art . The reconstruction efficiency is different for various R. 



1.4 Inclusive cross section of jet production 

We study the dependence of jet cross section on the jet-finding algorithm parameters 
E seed , R, E cut . Figure [5] presented the dependence of inclusive jet cross section on the 
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jet transverse energy E^ et for one, two, three and all-jet events (the left frame). The 
dependence of the cross section ratio on Ej, et is shown on the right frame. One see that 
events with one reconstructed jet dominate for Ej, et <15 GeV. Their contributions are 
noticeable for E^ et < 25 GeV. The dijet events dominate for E^ et > 15 GeV. The fraction 
of events with three reconstructed jets is noticeable over a range 20< E^ et <50 GeV. 
Maximum contribution of such events is about 15%. Behavior of jet multiplicity versus 
Ej, et is almost the same for another set of parameters (E cut = 5 GeV, E see d = 0.5 GeV, R = 
0.4). Maximum contribution of three-jet events is ~ 20%. 
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Figure 5: The dependence of the inclusive cross section (a) and the cross section ratio 
Jet 7v / Jet Aii (b) for all and iV-jet events (N = 1,2,3) on the jet transverse energy E^ et 
for PARI set of parameters in p — p interaction at ^=200 . (PARI : E^ et = 7 , E seed 
= 1.0,^ = 0.7). 




Figure 6: The dependence of the (Jet P AR2 — J ztpARi) I J etpAm ratio for two sets of 
parameters PARI PAR2 on the jet transverse energy E^ et for one-jet (a) and iV-jets 
(N Jet = 1, 2, 3) (b) events. (PARI: E#* = 7 , E seed = 1.0, R = 0.7). 

Figure M shows the dependence of the (Jet pari — JetpAm) / JetpAm ratio for two sets 
of parameters PARI PAR2 on the jet transverse energy E^f et . The left frame is for one-jet 
events. The right frame is for TV-jets (Nj et = 1,2,3). As seen from Figure [6] the ratio 
is practically constant above the threshold Ejf 1 > 15 GeV for dijet events. The absolute 
value of the obtained jet cross sections depend on the algorithm parameters. But the 
energy dependence of cross sections is practically the same for any sets of parameters 



5 



above E^ h = (2 — 3)E cut . Below this threshold the shape of cross sections depend very 
strongly on parameter choice. The value of the threshold for inclusive jet cross sections 
is found to be E$ h = (3 - 5)E cut . 

Thus, above the threshold E§, h the parameter choice results in the cross section normal- 
ization only. The value of cross section depends very strongly on the algorithm parameters 
below threshold. 



2 ^-Presentation of jet cross section 

In this section we would like to remind the basic ideas of z-scaling 0, Hj dealing with 
the investigation of the inclusive process 

Pi + P2 -> V + X. (4) 

Here the momenta and masses of the colliding and inclusive particles are denoted by 
Pi, P2,p and Mi, M 2 , mi, respectively. It is assumed that the process can be described in 
terms of the corresponding kinematic characteristics of the exclusive subprocess written 
in the symbolic form 

(xiMi) + {x 2 M 2 ) -»mi + {x 1 M 1 + x 2 M 2 + m 2 ). (5) 

The parameter m 2 is introduced to take into account the internal conservation laws (for 
isospin, baryon number and strangeness). The quantities x\ and x 2 are the fractions of 
incoming four-momenta Pi,P 2 of colliding objects. The energy of the parton subprocess 
is defined as 

s 1/2 = {x\M\ + 2x 1 x 2 (P 1 P 2 ) + x\Mlf 12 . (6) 

The elementary parton-parton collision is considered as a binary subprocess which satisfies 
the 4-momentum conservation law written in the following form 

{x x P x + x 2 P 2 - p) 2 = {x 1 M 1 + x 2 M 2 + m 2 ) 2 . (7) 

To determine the fractions X\ and x 2 the principle of minimum resolution is used. It 
states that the measure of the constituent interaction which connects kinematic and 
dynamic characteristics of the interaction 

n- 1 (x 1 ,x 2 ) = [(l-x 1 ) s '(l-x 2 ) 5 T\ (8) 

should be minimum one under the condition ([7]). Here 5i and S 2 are the factors relating 
the fractal structure of colliding objects [3j H]. 

The scaling variable z and scaling function ip{z) as suggested in (3J, S] are determined 
as follows _ 

§1 n-i if \ ns T-ij7 d3a fn\ 
ZqVL , ip(z) = J E—. (9) 



mp(s)Q ' ' p(s,rj)a in i dp 3 

Here m is the mass constant (nucleon mass), is the transverse energy of parton 

subprocess. 

The scaling variable z is expressed via the transverse kinetic energy of subprocess 
the multiplicity density of charge particles p(s) and Q(xi,x 2 ). The scaling function 
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i/j(z) is expressed via the colliding energy y/s, the average charged particle multiplicity 
density p(s,T]) = d < N > /drj, the inelastic cross section a in [, the inclusive cross section 
Ed 3 a/dp 3 , and the Jacobian J of transformation from the variables {p z ,Pr} to {i],z}. 
The function ip is interpreted as the probability density to produce a particle with the 
formation length z. Authors of z-scaling concept argue that the z-scaling reflects features 
of particle substructure, constituent interaction and particle formation such as locality, 
self-similarity and fractality over a wide scale range. 

The properties of the scaling function ip of jet production at ISR, SppS, and Tevatron 
were investigated in [31 H]. The energy and angular independence of ip and power behavior, 
if) ~ z _/3 , of jet and dijet production at high energy were established. It was found that 
the values of the slope parameter (3 in p — p and p — p collisions are different. The 
properties of the scaling function ip were used to predict the transverse energy spectrum 
of jet production in p — p and p — p collisions at RHIC, Tevatron and LHC energies. 

The STAR collaboration presented [6j the first results on inclusive cross sections of 
jet production in p — p collisions at RHIC. The transverse spectra are measured over the 
range = 200 GeV, 0.2 < \rj\ < 0.8, and 5 < E J T et < 50 GeV. 

In the present paper we analyze MC-results and RHIC data [6] in z-presentation. The 
comparison with predictions of z-scaling is given as well. 
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Figure 7: The inclusive cross section of jet production in p — p collisions at y/s = 200 in 
Et (a,b) and z (c) presentation. 

Figure [7] (a) shows Monte Carlo results on the invariant inclusive jet cross section for 
sets of algorithm parameters PARI and PAR2. Predictions of z-scaling are shown by the 
solid line. We found that the shape of MC jet cross sections for different sets of parameters 
coincides with the curve predicted by z-scaling over a range 25 < E^ et < 55 GeV. As seen 
from [7] (b) the STAR date are in agreement with MC results. The strong dependence of 
the cross section shape on the algorithm parameters with decreasing E^ et is observed for 
E^ et < 25 GeV. 

The comparison of the MC results, STAR data and predictions of z-scaling in z- 
presentation are shown in[7^c). We see that both MC simulation and STAR data are in a 
good agreement with z-scaling predictions for E^ et > 25 GeV (z > 180). The value of the 
slope parameter j3 is found to be (3 = 6.01 ±0.06 for Monte Carlo results at -^i = 200 GeV. 
It is compatible with f3 = 5.95 ± 0.21 obtained for the AFS data at -^/s = 38.8,45, and 
63 GeV [1]. Note that the shape of the scaling function ip(z) for p T < 10 GeV/c can not 
be described by the power law ip(z) ~ z _/3 . The error bars of the experimental data [6] 
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are large enough for precise test of asymptotic behavior of ^-scaling of jet production in 
p — p collisions at RHIC energies. 



3 Conclusions 

The Monte Carlo study of jet production in p — p collisions at y/s = 200 GeV using 
the PYTHIA generator was performed. The impact of the one algorithm parameters 
E cut , E see di R on the efficiency and characteristics of the reconstructed jets was investi- 
gated. It was established that the P^ art reconstruction efficiency is independent of E cut 
at R = 0.7. and depends on R for two leading jets in A-jet events. 

The probability of two-jet reconstruction was found to drop very fast if the transverse 
energy of hard process is less them the threshold p± < PJ^ m . It goes to saturation 
if p± > -Pf"*. These limits are controlled by the values of the algorithm parameters 
E seed , R, Ecut. Results of Monte Carlo study were compared with predictions made in 
the framework of z-scaling and experimental data obtained by the STAR collaboration at 
RHIC. The independence of the slope parameter f3 of the scaling function, ip(z) z~@, on 
the algorithm parameters over a range E^ et = 25 — 60 GeV was established. The strong 
dependence of the transverse energy spectra of jet production and the slope parameter 
on algorithm parameters were found for E^ et < 25 GeV. 

Verification of ^-scaling for jet production in p — p collisions at RHIC and compar- 
ison with data at Tevatron and LHC energies are of interest for understanding the jet 
phenomena. 
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